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Abstract 

In this paper we describe a dynamic data structure that answers one-dimensional stabbing- 
max queries in optimal 0(logn/ log log n) time. Our data structure uses linear space and sup- 
ports insertions and deletions in O(logn) and 0(logn/loglogn) amortized time respectively. 

We also describe a (9(n(logn/loglogn) d_1 ) space data structure that answers d-dimensional 
stabbing-max queries in 0((logn/loglogn) d ) time. Insertions and deletions are supported in 
0((logn/loglogn) £i loglogn) and 0((log n/ log log n) d ) amortized time respectively. 

1 Introduction 

In the stabbing-max problem, a set of rectangles is stored in a data structure, and each rectangle s 
is assigned a priority p(s). For a query point q, we must find the highest priority rectangle s that 
contains (or is stabbed by) q. In this paper we describe a dynamic data structure that answers 
stabbing-max queries on a set of one-dimensional rectangles (intervals) in optimal O (log nj log log n) 
time. We also show how this result can be extended to d > 1 dimensions. 

Previous Work. The stabbing-max problem has important applications in networking and 
geographic information systems. Solutions to some special cases of the stabbing-max problem play 
a crucial role in classification and routing of Internet packets; we refer to e.g., [9| lilt [18] for a small 
selection of the previous work and to |X0^ [T9] for more extensive surveys. Below we describe the 
previous works on the general case of the stabbing-max problem. 

The semi-dynamic data structure of Yang and Widom [24] maintains a set of one-dimensional 
intervals in linear space; stabbing-max queries and insertions are supported in O(logn) time. Agar- 
wal et al. [1] showed that stabbing-max queries on a set of one-dimensional intervals can be answered 
in 0(log 2 n) time; the data structure of Agarwal et al. [1] uses linear space and supports updates in 
0(log n) time. The linear space data structure of Kaplan et al. [15] supports one-dimensional queries 
and insertions in O(logn) time, but deletions take O(lognloglogn) time. In [15] , the authors also 
consider the stabbing-max problem for a nested set of intervals: for any two intervals si, S2 G S, 
either s\ C S2 or S2 C si or si n S2 = 0. Their data structure for a nested set of one-dimensional 
intervals uses O(n) space and supports both queries and updates in O(logra) time. Thorup [22] 
described a linear space data structure that supports very fast queries, but needs log^ 1 ) n time to 
perform updates. His data structure supports stabbing-max queries in 0(£) time and updates in 
0{n 1 ^) time for any parameter / = o(logn/ log log n). However the problem of constructing a data 
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structure with poly- logarithmic update time is not addressed in [22J. Agarwal et al. [2] described 
a data structure that uses linear space and supports both queries and updates in O(logn) time for 
an arbitrary set of one-dimensional intervals. The results presented in [15] and [2] are valid in the 
pointer machine model [20j . 

The one-dimensional data structures can be extended to d > 1 dimensions, so that space 
usage, query time, and update time increase by a factor of 0(log d_1 n). Thus the best previously 
known data structure [2] for ci-dimensional stabbing- max problem uses 0(n log d_1 n) space, answers 
queries in 0(log d n) time, and supports updates in 0(log d n) time. Kaplan et al. [15] showed that 
d-dimensional stabbing-max queries can be answered in O(logn) time for any constant d in the 
special case of nested rectangles. 

Our Result. The one-dimensional data structure described in [2] achieves optimal query time 
in the pointer machine model. In this paper we show that we can achieve sublogarithmic query time 
without increasing the update time in the word RAM model of computation. Our data structure 
supports deletions and insertions in O (log n/ log log n) and 0(log n) amortized time respectively. As 
follows from the lower bound of [3], any fully-dynamic data structure with poly-logarithmic update 
time needs f2(logn/loglogn) time to answer a stabbing-max quer}Q Thus our data structure 
achieves optimal query time and space usage. 

Our result can be also extended to d > 1 dimensions. We describe a data structure that uses 
0(n(logn/ loglogn) d_1 ) space and answers stabbing-max queries in (9((logn/loglogn) rf ) time; 
insertions and deletions are supported in 0((logn/loglogn) d loglogn) and 0((logn/ loglogn) d ) 
amortized time respectively. Moreover, our construction can be modified to support stabbing-sumn 




queries: we can count the number of one-dimensional intervals stabbed by a query point q in 
0(logn/ log log n) time. The stabbing-sum data structure also supports insertions and deletions in 
O(logn) and 0(logn/loglogn) amortized time. 

Overview. We start by describing a simple one-dimensional stabbing-max data structure in 
section [2] This data structure achieves the desired query and update times but needs uj{n 2 ) space. 
All intervals are stored in nodes of the base tree of height 0(logn/loglogn); the base tree is 
organized as a variant of the segment tree data structure. Intervals in a node u are stored in a 
variant of the van Emde Boas (VEB) data structure [8] . We can answer a stabbing-max query by 
traversing a leaf-to-root path in the base tree; the procedure spends O(l) time in each node on the 
path. 

In section [3j we show how all secondary data structures in the nodes of the base tree can be 
stored in O(nlogn) bits of space. The main idea of our method is the compact representation of 
intervals stored in each node. Similar compact representations were also used in data structures 
for range reporting queries [6] [16] and some other problems [5] . However the previous methods are 
too slow for our goal: we need O(loglogn) time to obtain the representation of an element e in a 
node u if the representation of e in a child of u is known. Therefore it would take O(logn) time 
to traverse a leaf-to-root path in the base tree. In this paper we present a new, improved compact 
storage scheme. Using our representation, we can traverse a path in the base tree and spend 0(1) 
time in each node. We believe that our method is of independent interest and can be also applied 
to other problems. 

1 In fact, the lower bound of [3] is valid even for existential stabbing queries: Is there an interval in the set S that 
is stabbed by a query point q? 

2 The stabbing-sum problem considered in this paper is to be distinguished from the more general stabbing-group 
problem, in which every interval is associated with a weight drawn from a group G. 
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Figure 1: Internal nodes of the base tree. Intervals a and 8 are stored in the set S(w), but are 
not stored in the set S(u). Interval 7 is stored in S(u), but 7 is not stored in S(w). Interval j3 
is stored in both S(u) and S(w). Moreover, a, (3, and 8 belong to the sets £23 (w), £33 (it;), and 
£33 (w) respectively. Intervals (3 and 7 belong to S%i(u) and £23(1*) respectively. 

The results for multi-dimensional stabbing-max queries and stabbing-sum queries are described 
in Theorems [2] and |3j an extensive description of these results is provided in Appendices C and D. 

2 A Data Structure with Optimal Query Time 

Base Tree. In this section we describe a data structure that answers stabbing-max queries in 
optimal time. Endpoints of all intervals from the set S are stored in the leaves of the base tree 
T. Every leaf of T contains @(log e n) endpoints. Every internal node, except of the root, has 
0(log £ n) children; the root has 0(log e n) children. Throughout this paper e denotes an arbitrarily 
small positive constant. The range rng(u) of a node u is an interval bounded by the minimal and 
maximal values stored in the leaf descendants of u. 

For a leaf node ui, the set S(ui) contains all intervals s, such that at least one endpoint of s 
belongs to U[. For an internal node it, the set S(u) contains all intervals s, such that rng{ui) C s 
for at least one child U{ of u but rng(u) <f_ s. See Fig. [I] for an example. Thus each interval is 
stored in 0(logre/loglogn) sets S(u). For every pair i < j, Sij(u) denotes the set of all intervals 
s £ S(u) such that rng(uf) C s for a child Uf of u if and only if i < f < j. For simplicity, we will 
sometimes not distinguish between intervals and their priorities. 

Secondary Data Structures. For each internal node u, we store a data structure D{u) 
described in the following Proposition. 

Proposition 1 Suppose that priorities of all intervals in S(u) are integers in the interval [l,p max \ 
for p m ax = 0(n). There exists a data structure D{u) that uses 0(p max ■ log 2e n) words of space 
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and supports the following queries: for any I < r, the predecessor of q in S[ r (u) can be found 
in O(loglogn) time and the maximum element in Si r (u) can be found in O(l) time. The data 
structure supports insertions in O(loglogra) time. If a pointer to an interval s € S(u) is given, s 
can be deleted in 0(1) amortized time. 

Proof: It suffices to store all intervals from Si r (u) in a VEB data structure Di r (u). Each Di r (u) 
uses 0(p max ) words of space and answers queries in 0(loglogp max ) = O(loglogn) time [8]. It is 
a folklore observation that we can modify the VEB data structure so that maximum queries are 
supported in constant time. □ 

We also store a data structure M(u) that contains the interval maxjj(u) with maximal priority 
among all intervals in Sij(u) for each i < j. Since each internal node has G(log £ n) children, M(u) 
contains 0(log 2£ n) elements. For any query index /, M(u) reports the largest element among all 
maxjj, i < f < j. In other words for any child Uf of u, M(u) can find the interval with the highest 
priority that covers the range of the /-th child of u. Using standard techniques, we can implement 
M(u) so that queries and updates are supported in 0(1) time. For completeness, we describe the 
data structure M(u) in Appendix A. 

Queries and Updates. Let tt denote the search path for the query point q in the base tree T. 
The path tt consists of the nodes Vq, v\, . . . , vr where vq is a leaf node and vr is the root node. Let 
s(vi) be the interval with the highest priority among all intervals that are stored in Uj<iS(vj) and 
are stabbed by q. The interval s(vq) can be found by examining all 0(log e ra) intervals stored in 
S(vo). Suppose that we reached a node Vi and the interval s(vi-\) is already known. If q stabs an 
interval s stored in S(vi), then rng(vi-\) C s. Therefore q stabs an interval s 6 S(vi) if and only if 
s is stored in some set Si r (vi) such that I < f < r and Vi-\ is the /-th child of Vi. Using the data 
structure M(fj), we can find in constant time the maximal interval s m , such that s m G Si r (vt) and 
I < / < t. Then we just set s(vi) = max(s m , s(vi-\)) and proceed in the next node Vi + \. The total 
query time is O (log nj log log n). 

When we insert an interval s, we identify 0(log n/ log log n) nodes Vi such that s is to be inserted 
into S(vi). For every such Vi, we proceed as follows. We identify I and r such that s belongs to 
Si r (vi). Using D(vi), we find the position of s in S[ r (vi), and insert s into S[ r (vi). If s is the 
maximal interval in Si r (vi), we delete the old interval maxi r (vi) from the data structure M(vi), set 
maxi r (vi) = s, and insert the new max( r (i)j) into M(vi). 

When an interval s is deleted, we also identify nodes Vi, such that s £ S(vi). For each Vj, we 
find the indices / and r, such that s S S[ r (vi). Using the procedure that will be described in the 
next section, we can find the position of s in Si r (vi). Then, s is deleted from the data structure 
D(vt). If s = m&xi r (vi), we remove maxi r (vi) from M(vi), find the maximum priority interval in 
Si r (vi), and insert it into M(v{). We will show in section [3] that positions of the deleted interval s 
in Si r (vi) for all nodes Vi can be found in (9(logn/ log log n) time. Since all other operations take 
O(l) time per node, the total time necessary for a deletion of an interval is 0(logn/loglogn). 

Unfortunately, the space usage of the data structure described in this section is very high: every 
VEB data structure Di r (u) needs 0(p ma _ x ) space, where p max is the highest possible interval priority. 
Even if p max = 0(n), all data structures D(u) use 0(n 2 log 2e n) space. In the next section we show 
how all data structures D(u), u GT, can be stored in 0(n log n) bits without increasing the query 
and update times. 
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3 Compact Representation 



The key idea of our compact representation is to store only interval identifiers in every node u of 
T. Our storage scheme enables us to spend O (log log n) bits for each identifier stored in a node u. 
Using the position of an interval s in a node u, we can obtain the position of s in the parent w 
of u. We can also compare priorities of two intervals stored in the same node by comparing their 
positions. These properties of our storage scheme enable us to traverse the search path for a point 
q and answer the query as described in section [2} 

Similar representations were also used in space-efficient data structures for orthogonal range 
reporting [6], [16] and orthogonal point location and line-segment intersection problems [5]. Storage 
schemes of |16t [5] also use O (log log n) bits for each interval stored in a node of the base tree. The 
main drawback of those methods is that we need O(loglogn) time to navigate between a node 
and its parent. Therefore, G(logn) time is necessary to traverse a leaf-to-root path and we need 
fi(logn) time to answer a query. In this section we describe a new method that enables us to 
navigate between nodes of the base tree and update the lists of identifiers in O(l) time per node. 
The main idea of our improvement is to maintain identifiers for a set S(u) D S(u) in every u G T ■ 
When an interval is inserted in S(u), we also add its identifier to S(u). But when an interval is 
deleted from S(u), its identifier is not removed from S(u). When the number of deleted interval 
identifiers in all S(u) exceeds the number of intervals in S(u), we re-build the base tree and all 
secondary structures (global re-build). 

Compact Lists. We start by defining a set S'(u) D S(u). If u is a leaf node, then S'{u) = S(u). 
If u is an internal node, then S'(u) = S(u) U (UjS"(ui)) for all children iij of u. An interval s belongs 
to S'(u) if at least one endpoint of s is stored in a leaf descendant of u. Hence, | U v S'(v)\ = 0(n) 
where the union is taken over all nodes v that are situated on the same level of the base tree T ■ 
Since the height of T is 0(logn/loglogn), the total number of intervals stored in all S'(u), u G 7~, 
is O (n log n / log log n) . 

Let S(u) be the set that contains all intervals from S'(u) that were inserted into S'(u) since the 
last global re-build. We will organize global re-builds in such way that at most one half of elements 
in all S(u) correspond to deleted intervals. Therefore the total number of intervals in \J U £tS(u) 
is 0{n log nj log log n). We will show below how we can store the represntations of sets S(u) in 
compact form, so that an element of S(u) uses O(loglogn) bits in average. Since S(u) C S(u), we 
can also use the same method to store all S(u) and D (u) in 0(n log n) bits. 

Sets S'(u) and S(u) are not stored explicitly in the data structure. Instead, we store a list 
Comp(u) that contains a compact representation for identifiers of intervals in S(u). Comp(u) is 
organized as follows. The set S(u) is sorted by interval priorities and divided into blocks. If 
> log 3 n/2, then each block of S(u) contains at least log 3 re/2 and at most 21og 3 n elements. 
Otherwise all e G S(u) belong to the same block. Each block B is assigned an integer block label 
lab(i?) according to the method of |14|. I23j . Labels of blocks are monotone with respect to order of 
blocks in Comp(n): The block B\ precedes B2 in Comp(u) if and only if lab(i^i) < lab(i?2). Besides 
that, all labels assigned to blocks of Comp(n) are bounded by a linear function of |Comp(n)|/log 3 n: 
for any block B in Comp(u), lab(S) = 0(|Comp(n)|/log 3 n). When a new block is inserted into a 
list, we may have to change the labels of 0(log 2 n) other blocks. 

For every block B, we store its block label as well as the pointers to the next and the previous 
blocks in the list Comp(u). For each interval s in a block B of S(u), the list Comp(u) contains the 
identifier of s in Comp(u). The identifier is simply the list of indices of children itj of u, such that 
s G Comp(uj). As follows from the description of the base tree and the sets S(u), we store at most 
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two child indices for every interval s in a block; hence, each identifier uses O(loglogn) bits. To 
simplify the description, we sometimes will not distinguish between an interval s and its identifier 
in a list Comp(u). 

We say that the position of an interval s in a list Comp(u) is known if the block B that contains 
s and the position of s in B are known. If we know positions of two intervals s\ and S2 in Comp(u), 
we can compare their priorities in 0(1) time. Suppose that s\ and S2 belong to blocks B\ and 
B2 respectively. Then s\ > S2 if lab(Si) > lab(2?2), and s\ < S2 if lab(i?i) < lab(.E?2). If 
lab(i?i) = lab(i?2), we can compare priorities of s\ and S2 by comparing their positions in the 
block Bi = B 2 . 

The rest of this section has the following structure. First, we describe auxiliary data structures 
that enable us to search in a block and navigate between nodes of the base tree. Each block B 
contains a poly-logarithmic number of elements and every identifier in B uses O(loglogn) bits. We 
can use this fact and implement block data structures, so that queries and updates are supported 
in O(l) time. If the position of some s in a list Comp(u) is known, we can find in constant time the 
positions of s in the parent of u. Next, we show how data structures D(u) and M(u), defined in 
section [2j are modified. Finally, we describe the search and update procedures. 

Block Data Structures. We store a data structure F{B) that supports rank and select queries 
in a block B of Comp(n): A query rank(/, i) returns the number of intervals that also belong to 
S(uf) among the first i elements of B. A query select(/, i) returns the smallest j, such that 
rank(/, j) = i; in other words, select(/, i) returns the position of the i-th interval in B that also 
belongs to S(uf). We can answer rank and select queries in a block in O(l) time. We can also 
count the number of elements in a block of Comp(u) that are stored in the /-th child of u, and 
determine for the r-th element of a block in which children of u it is stored. Implementation of 
F(B) will be described in Appendix A. 

For each block Bj £ Comp(u) and for any child Uf of u, we store a pointer to the largest block 
Bj before Bj that contains an element from Comp(?ij). These pointers are maintained with help 
of a data structure Pf(u) for each child Uf. We implement Pf(u) as an incremental split-find data 
structure |13j . Insertion of a new block label into Pf(u) takes O(l) amortized time; we can also 
find the block Bj for any block Bj £ Comp(u) in 0(1) worst-case time. Using the data structure 
F(Bj) for a block Bj, we can identify for any element e in a block Bj the largest e* < e, ej £ Bj, 
such that ej belongs to Comp(uy). Thus we can identify for any e £ Comp(u) the largest element 
ef £ Comp(u), such that e/ < e and e/ £ Comp(uj), in 0(1) time. 

Navigation in Nodes of the Base Tree. Finally, we store pointers to selected elements in 
each list Comp(u). Pointers enable us to navigate between nodes of the base tree: if the position of 
some e £ Comp(u) is known, we can find the position of e in Comp(u;) for the parent w of u and the 
position of e in Comp(uj) for any child Ui of u such that Comp(uj) contains e. 

We associate a stamp t(e) with each element stored in a block B; every stamp is a positive 
integer bounded by 0(|B|). A pointer to an element e in a block B consists of the block label 
lab(i?) and the stamp of e in B. When an element is inserted into a block, stamps of other 
elements in this block do not change. Therefore, when a new interval is inserted into a block B 
we do not have to update all pointers that point into B. Furthermore, we store a data structure 
H(B) for each block B. Using H(B), we can find the position of an element e in B if its stamp 
t(e) in B is known. If an element e must be inserted into B after an element e p , then we can assign 
a stamp t e to e and insert it into H{B) in 0(1) time. Implementation of H{B) is very similar to 
the implementation of F{B) and will be described in the full version. 
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Figure 2: Parent and child pointers in blocks of Comp(n) and Comp(-Uj). Intervals that are stored 
in Comp(uj) are depicted by circles; filled circles are intervals with pointers. Only relevant intervals 
and pointers are shown in Comp(-u). In Comp(uj), only parent pointers are shown. 

If e is the first element in the block B of Comp(-u) that belongs to Comp(uj), then we store the 
pointer from the copy of e in Comp(u) to the copy of e in Comp(uj). If an element e E Comp(u) is 
also stored in Comp(nj) and e is the first element in a block B' of Comp(iij), then there is a pointer 
from the copy of e in Comp(ii) to the copy of e in Comp(iij). Such pointers will be called child 
pointers. For any pointer from e E Comp(-u) to e G Comp(uj), we store a pointer from e G Comp(ui) 
to e € Comp(n). Such pointers will be called parent pointers. See Fig. [2] for an example. 

We can store each pointer in O(logra) bits. The total number of pointers in Comp(n) equals to 
the number of blocks in Comp(u) and Comp (it j) for all children Uj of u. Hence, all pointers and all 
block data structures use 0(n log n) bits. 

If we know the position of some interval s in Comp(t;), we can find the position of s in the parent 
w of v as follows. Suppose that an interval s is stored in the block B of Comp(u) and v is the f-th 
child of w. We find the last interval s' in B stored before s, such that there is a parent pointer 
from s'. Let m denote the number of elements between s' and s in B. Let B' be the block in 
Comp(u;) that contains s' . Using H(B'), we find the position m! of s' in the block B 1 of Comp(u)). 
Then the position of s in B' can be found by answering the query select (/, rank(/, m!) + m) to a 
data structure F(B'). Using a symmetric procedure, we can find the position of s in Comp(v) if its 
position in Comp(w) is known. 

Root and Leaves of the Base Tree. Elements of S(vr) are explicitly stored in the root 
node vr. That is, we store a table in the root node vr that enables us to find for any interval s 
the block B that contains the identifier of s in Comp(u#) and the stamp of s in B. Conversely, if 
the position of s in a block B of Comp(u/j) is known, we can find its stamp in B in O(l) time. If 
the block label and the stamp of an interval s are known, we can identify s in 0(1) time. In the 
same way, we explicitly store the elements of S(ui) for each leaf node ui. Moreover, we store all 
elements of S(vr) in a data structure R, so that the predecessor and the successor of any value x 
can be found in 0(logn/loglogn) time. 

Data Structures M(u) and D(u). Each set S(u) and data structure D(u) are also stored in 
compact form. D{u) consists of structures Di r (u) for every pair I < r. If a block B contains a label 
of an interval s G Si r (u), then we store the label of B in the VEB data structure Di r (u). The data 
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structure G{B) contains data about intervals in S(u) Pi B. For every s G G(B), we store indices 
l,r if s G Si r (u); if an element s £ B does not belong to 5 (i.e., s was already deleted), then we 
set Z = r = 0. For a query index /, G(B) returns in 0(1) time the highest priority interval s € B, 
such that s G Si r (u) and I < f < r. A data structure G(.B) uses 0(|J3| log log n) bits of space and 
supports updates in O(l) time. Implementation of G(B) is very similar to the implementation of 
F(B) and will be described in the full version of this paper. We store a data structure M{u) in 
every node u G T, such that Comp(u) consists of at least two blocks. For each pair I < r, the data 
structure M(u) stores the label of the block that contains the highest priority interval in S[ r (u). 
For a query /, M{u) returns the label of the block that contains the highest priority interval in 
Ui<f< r Si r (u). Such queries are supported in 0(1) time; a more detailed description of the data 
structure M{u) will be given in Appendix A. 

Search Procedure. We can easily modify the search procedure of section [2] for the case when 
only lists Comp(u) are stored in each node. Let Vi be a node on the path n = vq, . . . vr, where tt is 
the search path for the query point q. Let s(vi) denote the interval that has the highest priority 
among all intervals that belong to Uj<iS(vj) and are stabbed by q. We can examine all intervals 
in S(vq) and find s(vq) in 0(log e n) time. The position of s(vo) in Comp(t>o) can also be found in 
O(l) time. During the i-th step, i > 1, we find the position of s(i>i) in Comp(fj). An interval s in 
S(vi) is stabbed by q if and only if s G Sfr(uj), I < f < r, and is the f-th child of u$. Using 
M(vi), we can find the block label of the maximal interval s m among all intervals stored in Si r (vi) 
for I < f < r. Let B m be the block that contains s m . Using the data structure G(B m ), we can find 
the position of s m in B m . 

By definition, s(i)j) = max(s m , s(uj_i)). Although we have no access to s m and s(uj_i), we can 
compare their priorities by comparing positions of s m and s(«j_i) in Comp(wj). Since the position of 
s(vi-\) in Comp(?jj_i) is already known, we can find its position in Comp(vj) in O(l) time. We can 
also determine, whether s m precedes or follows s(uj_i) in Comp(wj) in O(l) time. Since a query to 
M(vi) also takes O(l) time, our search procedure spends constant time in each node vi for i > 1. 
When we know the position of s(vr) in Comp(v/j), we can find the interval s(vr) in O(l) time. The 
interval s(vr) is the highest priority interval in S that is stabbed by q. Hence, a query can be 
answered in 0(logn/ log log n) time. 

We describe how our data structure can be updated in section |4j Our result is summed up in 
the following Theorem. 

Theorem 1 There exists a linear space data structure that answers orthogonal stabbing-max 
queries in 0(logn/ log log n) time. This data structure supports insertions and deletions in O(logra) 
and 0(logn/loglogn) amortized time respectively. 

4 Updates in the Data Structure of Theorem [I] 

Suppose that an interval s is inserted into S. The insertion procedure consists of two parts. First, 
we insert s into lists Comp(-Uj) and update Comp(u,) for all relevant nodes V{ of T. Then, we insert 
s into data structures D{vj) and M[vj). 

Using the data structure R, we can identify the segment s'(vji) that precedes s in S(vr). Then, 
we find the position of s'(vr) in Comp(^). We also find the leaves v a and Vf, of T, in which the left 
and the right endpoints of s must be stored. 

The path ir = vr, . . . ,v a is traversed starting at the root node. Let s'(vi) be the segment that 
precedes s in Comp(wj) and let B{vi) be the block that contains s'{vi). In every node Vi, we insert 
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the identifier for s after s'(uj). We also update data structures F(B(vi)) and H(B(vi)) for the 
block B(vi) that contains s. If the number of intervals in B(vi) equals 21og 3 ra, we split the block 
B(vi) into two blocks B\{vi) and B%(yi) of equal size. We assign a new label to .62(^1) and update 
the labels for some blocks of Comp(t;j) in 0(log 2 n) time. We also may have to update 0(log £ n) 
split-find data structures Pf(vi). Besides that, 0(1) child pointers in the parent of V{ and 0(log £ n) 
parent pointers in the children of V{ may also be updated. We split the block after 0(log 3 n) 
insertions; hence, an amortized cost of splitting a block is 0(1). When s is inserted into Comp(i»j), 
we find the largest element s'(vi-i) < s, such that s'(vi-i) is also stored in Comp(uj_i). Then, we 
find the position of s'(uj_i) in Comp(wj_i) and proceed in the node 

Let v c be the lowest common ancestor of v a and Vb, and let 7r& be the path from v c to Vb- We 
also insert s into Comp(fj) for all Vj G 7T&; nodes Vj G 7Tfo are processed in the same way as nodes 

During the second stage, we update data structures D(vi) and M(vi) for v% G 7r a U 717,. For any 
such we already know the block B G Comp(vj) that contains s and the position of s in Hence, 
a data structure G{B) can be updated in 0(1) time. If s is the only interval in B that belongs 
to Si r (vi) for some I, r, we insert lab(S) into Di r (vi). If lab(£?) is the greatest label in D[ r (vi), 
we also update M(vi). An insertion into a data structure Di r (vi), Vi £ n a U -Kb, takes O(loglogra) 
time. Updates of all other structures in Vi take O(l) time. Hence, an insertion takes O(logn) time 
in total. 

When an interval is deleted, we identify its position in every relevant node Vi. This can be 
done in the same way as during the first stage of the insertion procedure. Then, we update the 
data structures D(vi), M(vi), and G(B(vi)) if necessary. Since a block label can be removed from 
Dirivi) in O(l) time, the total time necessary to delete an interval is 0(logn/loglogn). 

Re-Balancing of the Tree Nodes. It remains to show how the tree can be rebalanced after 
updates, so that the height of T remains O (log nj log log n). We implement the base tree T as a 
weight-balanced B-tree |5j with branching parameter <fi and leaf parameter <j> for <j) = log 6 n. We 
assume that the constant e < 1/8. 

When a segment s is deleted, we simply mark it as deleted. After n/2 deletions, we re- build 
the base tree and all data structures stored in the nodes. This can be done in 0(n log nj log log n) 
time. 

Now we show how insertions can be handled. We denote by the weight of u the total number 
of elements in the leaf descendants of u. The weight n u of u also equals to the number of segment 
identifiers in Comp(u). Our weight-balanced B-tree is organized in such way that n u = 0(^/ +1 ), 
where (j> = log £ n and i is the level of a node u. A node is split after Q(<p i+1 ) insertions; when a node 
u is split into u' and u", the ranges of the other nodes in the base tree do not change; see [1] for a 
detailed description of node splitting. We will show that all relevant data structures can be re-built 
in 0(n u ) time. Hence, the amortized cost of splitting a node is O(l). When a segment is inserted, 
it is inserted into 0(logn/loglogra) nodes of the base tree. Hence, the total amortized cost of all 
splitting operations caused by inserting a segment s into our data structure is O (log n j log log n). 

It remains to show how secondary data structures are updated when a node u is split. Let w 
be the parent of u. The total number of segments in S(u') and S(u") does not exceed |i?(it)| = 
0(n u ). Hence, we can construct data structures D(u'), D(u") and lists Comp(i/), Comp(u") with 
all block data structures in 0(|5(u)|) = 0{n u ) time. We can find the positions of all elements 
e G Comp(n / ) U Comp(u // ) in Comp(u)), update their identifiers in Comp(u;), and update all auxiliary 
data structures in 0(n u ) time. 
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Some of the intervals stored in S(u) can be moved from S(u) to S(w): if an interval s S S(u) 
does not cover rng(u) but covers rng(u') or rng(u"), then s must be stored in S(w) after splitting. 
See Fig. [3]for an example. The total number of elements in Comp(u;) is ®(<ft e+2 ) = @(n u ■ <j))\ hence, 
the total number of blocks in Comp(ii;) is o(n u / log 3 n) . Identifiers of all segments in S(u) are already 
stored in Comp(u> ). We can update all block data structures in 0(1) time per segment. Every update 
of the data structure D{w) takes O(loglogn) time. But the total number of insertions into D(w) 
does not exceed the number of blocks in D{w). Therefore the total time needed to update D{w) is 
<3(n u log £ ~ 3 n) = o(n u ). 

Suppose that u was the fc-th child of w. Some segments stored in Sik(w) for I < k can be 
moved to 5^fc + i) (w), and some segments in Sk r (w) for r > k can be moved to S^ + iy(w). Since 
the total number of blocks in Comp(u;) is 0(n u / log 3 n), all updates of data structures Dij(w) take 
0((n u / log 3 n) log log n) = o(n u ) time. At most n u segments were stored in Sik(w) and Sk r (w) 
before the split of u. Hence, we can update all block data structures in 0(n u ) time. 

Finally, we must change the identifiers of segments stored in Comp(u;) and data structures 
G(B). Recall that an identifier indicates in which children of w a segment s is stored. Suppose 
that a segment s had an identifier r > k in Comp(u7). Then its identifier will be changed to r + 1 
after the split of u. The total number of segments with incremented identifiers does not exceed 
n w = 0(n u ■ (p). However, each identifier is stored in O(loglogn) bits. We can use this fact, and 
increment the values of -^/log n identifiers in O(l) time using the following look-up table T. There 
arc 

0((l ogn )VI^) different 

sequences of y/\ogn identifiers. For every such sequence a and any 
j = 0(log £ n), T[a][j] = a' . Here a' is the sequence that we obtain if every f > j in the sequence 
a is replaced with j' + 1. Thus the list Comp(w) can be updated in 0(n w /\/logn) = 0(n u ) time. 
Using the same approach, we can also update data structures F(B) and G(B) for each block B in 
0(\B\/ \J\og n) time. Hence, the total time necessary to update all data structures in w because 
some identifiers must be incremented is 0{n u ). 

Since all data structures can be updated in 0(n u ) time when a node u is split, the total amortized 
cost of an insertion is O(logn). 
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5 Multi-Dimensional Stabbing-Max Queries 



Our data structure can be extended to the case of d-dimensional stabbing- max queries for d > 1. 
In this section we prove the following theorem. 

Theorem 2 There exists a data structure that uses 0(n(logn/loglogn) d_1 ) space and answers 
d- dimensional stabbing-max queries in 0((logn/loglogn) d ) time. Insertions and deletions are 
supported in 0((logn/loglogn) d loglogn) and 0((logn/loglogn) d ) amortized time respectively. 

Let (d, g)-stabbing problem denote the (d+g)-dimensional stabbing-max problem for the case when 
the first d coordinates can assume arbitrary values and the last g coordinates are bounded by log p n, 
for p < 8 fg +( n ■ First, we show that the data structure of Theorem 1 can be modified to support 
(1, g)-stabbing queries for any constant g. Then, we will show how arbitrary d-dimensional queries 
can be answered. Throughout this section proji(s) denotes the projection of a rectangle s on the 
i-th coordinate. We denote by bi(s) and ei{s) the i-th coordinates of endpoints in a rectangle s, so 
that proji(s) = [h(s), ej(s)]. 

A Data Structure for the (1, g)-Stabbing Problem. The solution of the (1, (7)-stabbing 
problem is very similar to our solution of the one-dimensional problem. We construct the base tree 
T on the ^-coordinates of all rectangles. T is defined as in section [2j but we assume that e < p. 
Sets S(u), Si r (u), and S(u) are defined as in sections [2j[3j We define Si r [j\, . . . , j g , . . . ,j2g](u) as the 
set of rectangles s in Si r (u) such that b 2 {s) = j%, . . ., b g+1 (s) = j g , e 2 (s) = j g+1 , . . ., e g+1 (s) = j 2g . 

A rectangle identifier for a rectangle s stored in a list Comp(u) consists of one or two tuples 
with 2(7 + 1 components. The first component of each tuple is an index j, such that s belongs 
to the j-th child of u; remaining 2g components are the last g coordinates of the rectangle end- 
points. The data structure M(u) contains the maximum priority rectangle in Si r [ji, . . . , j2g](u) 
for each l,r,j\, . . . ,j 2 g- M{u) can find for any (f,h\,..., h g ) the highest priority rectangle in all 
S[ r [ji, . . . , j2 g ](u), such that I < f < r, and ji < hi < j g+ i for 1 < i < g. We can implement 
M{u) as for the one-dimensional data structure. Data structure D(u) consists of VEB structures 
Di r [ji, . . . , j2 g ](u). If a block B of Comp(u) contains a rectangle from Si r \ji, . . . ,j2g](u), then lab(i?) 
is stored in Ar[ji, • • • ,32g]{u). 

We can answer a query q = (q x , q±, . . . , q g ) by traversing the search path n = Vo, . . . , Vr for q x 
in the base tree T ■ As in Theorem [TJ we start at the leaf node v o and find the maximal rectangle 
stored in S(vi) that is stabbed by the query point q. The search procedure in nodes v\, . . . ,vr is 
organized in the same way as the search procedure in section [3j The rectangle s(vi) is defined as in 
section[3j If Vi-\ is the /-th child of Vi, we answer the query (/, q\, . . . , q g ) using the data structure 
M{vi). Then we visit the block B m returned by M(vi) and find the last rectangle s m in B m with 
an identifier (/, b 2 (s m ), b g+ i(s m ), e 2 (s m ), . . . , e ff+ i(s m )) such that b i+ i{s m ) < qt < e i+ i(s m ) for 
1 < i < 9- Finally we compare s m with the rectangle s(fi-i) by comparing their positions in 
Comp(uj). Thus we can find the position of s(vi) in Comp(vj). When we reach the root vr, s(vr) is 
the highest priority segment that is stabbed by q. 

Updates of the list Comp(n) and all auxiliary data structures are implemented as in section |4j 

Lemma 1 There exists a data structure that answers (1, g) -stabbing queries in 0(logn/ log log n) 
time and uses 0(n) space. Insertions and deletions are supported in O(logn) and 0(\ogn/ log log n) 
amortized time respectively. 

A Data Structure for the (d, g)-Stabbing Problem. The result for (1, g)-stabbing can be 
extended to (d, ^-stabbing queries using the following Lemma. 
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Lemma 2 Suppose that there is a 0(n(logn/ loglogn) d ~ 2 ) space data structure D\ that answers 
(d — 1, <7 + 1)- stabbing queries in 0((logn/loglogn) ci ~ 1 ) time; D\ supports insertions and deletions 
in 0((logn/ loglogn)^ 1 loglogn) and 0{ (log n/ log log n) d ~ l ) amortized time respectively. 
Then there exists a 0(n(logn/loglogn) d_1 ) space data structure D2 that answers (d, g)- stabbing 
queries in 0((logn/loglogn) d ) time; D2 supports insertions and deletions in amortized time 
0((logn/ log log n) loglogn) and 0((logn/ log log n) d ) respectively. 

The main idea is to construct the base tree T on the d-th coordinates of rectangles and store 
a data structure for (d — l,g + l)-stabbing queries in each tree node. The tree is organized as in 
section [2] and the first part of this section. Let q^ denote the d-th coordinate of a point. Leaves 
contain d-ih coordinates of rectangle endpoints. A rectangle s is stored in a set S(u) for a leaf u if 
projd(s) n rng(u) 7^ and rng(u) <f_ projd(s). A rectangle s is stored in a set S(u) for an internal 
node u if rng(u) (jL projd(s) and rng(ui) C projd(s) for at least one child Ui of u. Sij(u) is the set 
of all intervals s 6 S(u) such that rng(uf) C projd(s) for a child uj of u if and only if i < f < j. 
We store a data structure T>(u) for (d — l,g + l)-stabbing queries in each internal node u. For a 
rectangle s £ S[ r (u), V(u) contains a rectangle s' such that proj g (s) = proj g (s') for g 7^ d and 
projd(s') = [l,r]. In other words, we replace the d-th coordinates of s's endpoints with I and r. 

A query q is answered by traversing the search path for q. For a leaf node Uo, we examine all 
rectangles in S{v$) and find the highest priority rectangle in S{vq). In an internal node Vi, i > 1, 
we answer the query g('Ui) = ((ft, ... , g^-i, /, qd+i-, ■ ■ ■ > fti+g) using the data structure V{vi). The 
query is obtained from q by replacing the d-th coordinate with an index /, such that is 
the /-th child of Uj. When a node v% is visited, our procedure finds the highest priority rectangle 
s(vi) that is stored in S(vi) and is stabbed by q. All rectangles that are stabbed by q are stored in 
some S(vi). Hence, the highest priority rectangle stabbed by q is the maximum rectangle among 
all s(vi). Since our procedure spends O((logn/loglogn) o!_1 ) time in each node, the total query 
time is 0((logn/ loglogn) rf ). 

When a rectangle s is inserted or deleted, we update 0(logn/ loglogn) data structures T>(u) 
in which s is stored. Hence, deletions and insertions are supported in 0((logn/ loglogn)^) and 
((log n/ log log n) d loglogn) time respectively. We can re-balance the base tree using a procedure 
that is similar to the procedure described in section [4} 

Proof of Theorem [2[ Theorem [2] follows easily from Lemmas [l] and [2j By Lemma [TJ there 
exists a linear space data structure that answers (l,d — l)-stabbing queries in 0(logn/ loglogn) 
time. We obtain the main result for ci-dimensional stabbing-max queries by applying the result of 
Lemma[2]to the data structure for (1, d — l)-stabbing queries. 

6 Stabbing-Sum Queries 

We can also modify our data structure so that it supports stabbing-sum queries: count the number 
of intervals stabbed by a query point q. 

Theorem 3 There exists a linear space data structure that answers orthogonal stabbing-sum 
queries in 0(logn/ loglogn) time. This data structure supports insertions and deletions in O(logn) 
and 0(logn/ loglogn) amortized time respectively. 

The only difference with the proof of Theorem [T] is that we store data structures for counting 
intervals instead of M(u). We maintain a data structure X(u) in each internal node u. For a query 
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index /, X{u) reports in O(l) time the total number of intervals in U;< j< r 5; r (u). In every leaf node 
ui, we maintain a data structure Z(ui) that supports stabbing sum queries on S(ui) and updates 
in O (log log n) time. 

As above, let it = vq . . . vr denote the search path for a query point q. In every node Vi E n, 
i > 1, we count the number n(v{) of intervals in Ui<f< r Si r (u) using X{vi); the index / is chosen 
so that Vi-i is the f-th child of V{. We also compute the number n(vo) of intervals that belong to 
S(vo) and are stabbed by q using Z{vq). An interval s is stabbed by q either if s is stored in S(vi), 
i > 1, and rng(v{-i) C s or if s is stored in S(vq) and s is stabbed by q. Hence, q stabs YIvg-k n ( v i) 
intervals. Each n{vi), i > 1, can be found in O(l) time and u{vq) can be found in O(loglogn) time. 
Hence, a query can be answered in 0(logn/loglogn) time. 

Now we turn to the description of the data structure X(u). Following the idea of pTj, we store 
information about the recent updates in a word B; the array A reflects the state of the structure 
before recent updates. A is a static array that is re-built after log 2£ n updates of X(u). When 
we construct A, we set A[f] = J2i<f< r \^lr( u )\- The word B contains one integer value m(l,r) for 
each pair / < r (m(Z, r) can also be negative, but the absolute value of each m(l,r) is bounded by 
0(log e n)). When a segment is inserted into (deleted from) Si r (u), we increment (decrement) the 
value of m(l,r) by 1. We can find X^K/<r m (^ r ) ™ time using a look-up table. After log e n 
updates we rebuild the array A and set all m(/, r) = 0. The amortized cost of rebuilding A is 0(1). 

The total number of segments in U/<j< r S , ; r (ti) equals to A[f] + J2i<f< r m (^ r )- Hence, a 
query to X{u) is answered in 0(1) time. We implement A in such way that each entry A[i] 
uses 0(log|5(u)|) bits. Thus each data structure X(u) uses 0(log £ nlog |5(w)|) bits and all X(u), 
u G T, use 0(n log n) bits in total. 

It remains to describe the data structure Z(u). Let E(u) be the set that contains endpoints 
of all intervals in S(u). Since u is a leaf node, S(u) contains 0(log 2e n) elements. Hence, it takes 
O(loglogra) time to find the rank r(q) of q in E(u). For each 1 < % < \E\, C[i] equals to the number 
of intervals that are stabbed by a point q with rank r(q) = i. The array C enables us to count 
intervals stabbed by q in 0(1) time if the rank of q in E(u) is known. Since \E\ = 0(log e n) and 
C[i] = 0(log e n), the array C uses 0(log £ nloglogn) bits of memory. When a set S(u) is updated, 
we can update C in 0(1) time using a look-up table. Thus Z(u) uses linear space and supports 
both queries and updates in O (log log n) time. 
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Appendix A. Auxiliary Data Structures 



Using bitwise operations and table look-ups, we can implement the data structure M(u) and the 
block data structures so that queries are supported in constant time. 

Data Structure M(u). Let A4(u) be the set of all interval priorities stored in M(u). Recall 
that Ai(u) contains an element maxjj for each set Sij(u), where maxjj is the highest priority interval 
(or its block label) stored in S{j(u). For simplicity, we assume that maxjj = — oo if Sij = 0. Let 
the rank of x in M.(v) be defined as rank(:r, Ai(u)) = |{ e 6 M.{u) \ e < x}\. Let Ai'(u) be the set 
in which every element of Ai(u) is replaced with its rank in Ai(u). That is, for each maxjj stored 
in M(u) we store max^ = rank(maxjj, M.(u)) in the set M.'{u). 

Each of ©(log 1 / 4 n) elements in A4.'{u) is bounded by Oilog 1 ^ n). Hence, we can store M.'(u) in 
one bit sequence W; W consists of 0(log 1 / 4 nloglogn) bits and fits into a machine word. We can 
store a table Tbli with an entry for each possible value of W and for each /. The entry TblfVF^/] 
contains the pair <l,r>, such that maxj r is the highest value in the set { max^ | i < f < j }. 

Obviously, max„(, < max^ if and only if max' ab < max'^. Therefore a query / can be answered 
by looking up the value </, r >= Tbl[IU][,/] for W = M'{u) and returning the element max;,.. 

When M(u) is updated, the value of some maxjj is changed. We store all elements of Ai(u) in 
an atomic heap p2J [21] that supports predecessor queries and updates in O(l) time. Hence, the 
new rank of maxjj in A4(u) can be found in O(l) time. If the rank of maxjj has changed, the ranks 
of all other elements in Ai(u) also change. Fortunately, max^ can assume only 0(log 1/,4 n) values. 
There are ©(log 1 / 4 ?!) elements max^ in a set Ai'(u), and there are 

2 0(io g / n) = Q ( n ) different 

sets Ai'(u). Besides that, each W = Ai'(u) fits into one word. Hence, we can update the set 
A4'(u) using a look-up in a table Tbli. Suppose that the rank of the /-th element in a set A4'(u) 
is changed to r. Then the new set Ai'(u) is stored in an entry Tbli[H^][/][r] of Tbli; here W is the 
bit sequence that corresponds to the old set Ai'(u). Thus an update of M(u) can be implemented 
in O(l) time. 

We need only one instance of tables Tbl and Tbli for all nodes u of the base tree. Both tables 
use o(n) space and can be initialized in o(n) time. 

Data Structure F(B). The data structure F(B) for a block B £ Comp(tt) is implemented as 
a B-tree TT>. Every leaf of TT> contains &(p) identifiers and each internal node has degree Q(p) 
for p = \/Togn. Recall that each identifier consists of at most two indices i\, %2 such that the 
corresponding interval is stored in the «i-th and Z2-th children of the node u in the base tree. All 
identifiers stored in a leaf node can be packed into O(ploglogn) bits. In every internal node v of 
Tp, we store the bit sequence W(v). For every child vi and for every possible value of j, W{v) 
contains the total number of indices j in the i-th child v% of v\ W{v) consists of 0(/jloglogn) bits. 
Using a look-up table, we can count for any W{y) and for any i the total number of elements stored 
in the children ui, . . . , v%^\ of v. For any W{y) and i, we can also count the total number of indices 
j in the children u\, . . . , i/j_i of v. For any k < 2 log 3 n and any W(v), j, we can find the largest i 
such that the total number of indices j in v\ . . . , V{ does not exceed k. Look-up tables that support 
such queries use o(n) space and can be initialized in o(n) time. We need only one instance of each 
table for all blocks B and all F(B). 

All queries and updates of a data structure F(B) can be processed by traversing a path in TT>. 
Using the above described look-up-tables, we spend only constant time in each node of TT>; hence, 
queries and updates are supported in 0(1) time. Details will be given in the full version of this 
paper. Data structures G(B) and H(B) can be implemented in a similar way. 
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